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Abstract 

Using recent values of the QCD (non-) perturbative parameters given in Table[T[ we reconsider the extraction of f B and the on-shell 
mass Mt, from HQET Laplace spectral sum rules known to N2LO PT series and including dimension 7 condensates in the OPE. 
We especially study the convergence of the PT series, the effects on "different spectral sum rules data" of the continuum threshold 
and subtraction point varied in a larger range than in the existing literature and include in the error an estimate of the N3LO PT 
series based on a geometric growth of the PT series. We obtain the Renormalization Group Invariant (RGI) universal coupling : 

^ fg = 0.407(28) GeV 3/2 in the static limit M b -> °o and the physical decay constant including l/M b corrections: f^ et = 197(13) 
MeV. Using the ratio of sum rules, we obtain the on-shell b quark mass Mb = 4871(30) MeV to order a 2 s from which we deduce 
the running mass rhb(Jhb) — 4234(56) MeV. The previous results are in good agreement with the ones from QCD spectral sum rules 
(QSSR) in full QCD to the same order f B = 206 ( 7 ) MeV and m~b(™b) qcd = 4236(69) MeV from the same channel. Using 

^\ in addition the value m,b(mb) = 4177(11) MeV from the T sum rules to order a] [2], we deduce the final averaged estimate from 
QSSR: f B ssr = 204(6) MeV and m b {m b ) qssr = 4181(11) MeV. 

Keywords: QCD spectral sum rules, meson decay constants, heavy quark masses, heavy quark effective theory. 
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The (pseudo)scalar meson decay constants fp are of prime 
interests for understanding the realizations of chiral symmetry 
in QCD. In addition to the well-known values of f„ = (130.4(2) 
MeV and f K = 156.1(9) MeV which control the light flavour 
chiral symmetries, it is also desirable to extract the ones of 
the heavy-light charm and bottom quark systems with high- 
accuracy. These decay constants are normalized through the 
matrix element: 



(o\jy X )\p) 



f P M 2 p . 



£Nj where: 



JL(x) = (m q + M Q )q(iy 5 )Q . 
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k> is the local heavy-light pseudoscalar current; q = d,s; Q = 
;_i c,b; P = D( s ), B (J )and where f P is related to the leptonic width: 
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where mi is the lepton mass and \Vq 9 \ the CKM mixing angle. 
In a recent analysis [1|, we have revised the extraction of these 
heavy-light decay constants in full QCD 0] using QCD spec- 
tral sum rules [4—9]. Here, we pursue the analysis by revisiting 
the determination of from HQET spectral sum rules. In so 
doing, we shall explicitly analyze the influence on the results 
of the subtraction point fi, of the continuum threshold t c . We 



shall also use (besides recent determinations of the QCD input 
parameters) the new precise value of nib from the Y sum rules 
O. In addition, we shall re-extract the meson-quark mass dif- 
ference using HQET sum rules from which we shall deduce the 
running fo-quark mass. 



2. HQET preliminaries 

HQET spectral sum rules have been initially used by Shuryak 
ifm using a non-relativistic Q version of the NS VZ lfl2l sum 
rules in the large Mj limit. Shuryak's sum rule has been ap- 
plied later on in HQET [13] by several authors lfT4H20"ll . The 
most important input in the analysis of HQET sum rule is local 
heavy-light quark axial-vector current of the full QCD theory 
which can be expressed as an OPE of the HQET operators O n 
in the inverse of the heavy quark mass: 



J M Ax,Mb) = 



C b {^,a s iti)\F A { X ,M b = oo) 

+ h Cn \T' aM ) K ' (4) 



= qy^y 5 ^, is the quark current in HQET built from 
a light antiquark field u and a properly normalized heavy quark 
field h v lfT3l . Cb, n are Wilson coefficients and Mb is the on-shell 
/7-quark mass. Using a non-covariant normalization of hadronic 
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states EH . one can define an universal coupling in the static 
limit: 



<0|^|P(v)> = — f m ^ 
V2 



(5) 



The coefficient function Cb (M^/fi, a s (l-i)) is obtained by requir- 
ing that HQET reproduces the full QCD theory at jj. = Mq. It 
has been obtained to order a s in IfTTll and to order a 2 in Il22l . It 
reads in the MS -scheme: 
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with : a s = a s jn. The HQET current acquires anomalous 
dimension, which reads to 0(a 2 s ) lf]~3l l23l l24l (in our normal- 
izations): 



y = 7i a s +j 2 a s + 

127 7tt 2 5 
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Therefore, the universal coupling scales as: 
- R h (M b ) - 

fstatW = n , : fsfat(M b ) , 



Rb(M) 



with: 



R b (fi) = (a s (jj.)) 
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where the two first coefficients of the /3 functions are: 



(7) 
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The universal coupling is connected to the physical decay con- 
stant as: 



f B y/M B = Cb(Mb)f sla ,(M h ) + 0(1 1 M b ) 



(11) 



It is also convenient to introduce the universal Renormalization 
Group Invariant (RGI) current and the associated coupling: 



P = Rbi^PiM) , f B = Rb(M)fstat(M) 
which we shall estimate in the following. 



(12) 



3. HQET spectral sum rules for f B in the static limit 

We shall be concerned with the universal 2-point-function|2] 

tl(q 2 = -Q 2 ) = i j d 4 x e ic ' x (0\J(x)jH0)\0) (13) 



2 In HQET Lorentz structure is unimportant and it is only the parity which 
counts such that we shall omit it in the following. 



for determining the coupling f B using QCD spectral sum rules 
(QSSR). Like in the case of the full theory, we can use either 
the Laplace (LSR) |H|25]|26): 



lim 

Q 2 , n — > oo 
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or the Q 2 = Moments sum rules (MSR) 0: 
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However, the use of the Q 2 = 0-moment sum rules is rather 
delicate as they do not have a proper infinite heavy quark mass 
limit. We shall not consider these sum rules here For the 
present analysis, it is convenient to introduce respectively the 
soft scale, the meson-quark mass-difference and the HQET 
Laplace sum rule variable: 



(q 2 -M 2 h ) 
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co= ' " , AM = — j— — - , T H = TM b , (16) 
M b M b ) 

where Mb is the on-shell quark mass and t is the usual LSR 
variable used in the full QCD theory and has the dimension of 
GeV~ 2 . As usual, we parametrize the spectral function using 
the Minimal Duality Ansatz (MDA): 

'■ " s ' • Ml) + "QCD cont."#(? - t c ), (17) 



Imn(f) 



where the accuracy for the sum rule approaches has been explic- 
itly tested from heavy quarkonia systems in |1 j. The perturba- 
tive (PT) expression of the spectral function has been evaluated 
to order 0(a s ) (NLO) in H3 and to 0(a 2 ) (N2LO) in ET]. It 
reads: 



Imf[p7-(w) = 
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with: La, ee \n(/j 2 /cj 2 ). We estimate the 0(a 3 s ) (N3LO) by as- 
suming the geometric growth of the PT series ||28l as a dual to 
the effect of a 1/q 2 term 11291 [301 which parametrizes the UV 
renormalon contributions. 

NP corrections up to dimension 7 condensates has been ob- 
tained in lfl6l [3TI . Including all the previous corrections the 
sum rule read for Mb — * oo: 



x Jo 



do e-" TH ImllpT (u) + NP 
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3 Some attempts to use Q 2 = 0-moment sum rules can be found in 1201 
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where: a> c = (t c - M 2 )IM h and 
NP(p) = -(uu){n) 
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4. The QCD input parameters 

The QCD parameters which shall appear in the following 
analysis will be the on-shell bottom quark mass Mj (we shall 
neglect the light quark masses q = u here and in the follow- 
ing), the light quark condensate (qq), the gluon condensates 
(8 2 G 2 ) = (g 2 G; v G>Z v } and <g 3 G 3 > = (g 3 f abc G; v G h vp G^), the 
mixed condensate (qgcrGq) = (qg(T f ' v (A a /2)G^ [y q) = M^iqq) 
and the four-quark condensate p(qq) 2 , where p 2 indicates 
the deviation from the four-quark vacuum saturation. Their val- 
ues are given in Table [T] We shall work with the running light 
quark parameters known to order a 3 ll6l 171 1321 : 



(qgcrGq)(jj) 



rh q , Q (-/3 ias r m x C(a s ) 

-A^(-j8i«,) 2/A /C(a s ) 

-Ml^i-^a^lda,). 
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fhq Q is the RGI quark mass, p q is spontaneous RGI light quark 
condensate ll33l . The QCD correction factor C(a s ) in the previ- 
ous expressions is numerically: 



C(a s ) = 1 + 0.895 la s + 1.3715a, + ... : n f 
= 1 + 1.1755a, + 1.5008a;; + ... : n f 

which shows a good convergence. We shall use: 

a s (M T ) = 0.325(8) => a s (M z ) = 0.1192(10) 



(22) 



(23) 



from r-decays 04l l35l . which agree perfectly with the world 
average 2012 J36][37j: 



a s {Mz) = 0.1184(7) . 



(24) 



The value of the running (qq) condensate is deduced from the 
value of (m u +m d )(2) = (7.9+0.6) MeV obtained in 11381139 1 and 
the well-known GMOR relation: (m u + m^Xuu + dd) = -m 2 f 2 . 
We shall use the value of the RGI spontaneous mass to order a s 
for consistency with the known a s correction in the OPE: 



/}„ = 251(6) MeV. 



(25) 



The values of the running MS mass m/,(»t/,) recently obtained in 
Ref. (2) from bottomium sum rules, will also be used[^] Using 
the relation between the running mt,(mt,) - 4177(1 1) MeV from 
the Y-systems [2| and the on-shell (pole) Mf, masses (see e.g. 
017] E7J, one can deduce to order a 2 : 



M h = 4804(50)„, ^ a s (M b ) = 0.2326(22) 



(26) 



where the error is mainly due to the one of a s . This large er- 
ror has to be contrasted with the precise value of the running 



mass, and can be an obstacle for a precise determination of f B 
from HQET at a given y.. On can see in Section 8 that a direct 
extraction of the on-shell mass from the HQET at the same a 2 
order leads to about the same value and error which is an (a 
posteriori) self-consistency check of the value and error used in 
Eq. (26 1 for the analysis. We are aware that the inclusion of 



the known a 3 -correction and an estimate of the PT higher or- 
der terms using a geometric growth of the PT coefficients a la 
Ref. J28]|29l increase the value of M b by about (100 ~ 200) 
MeV, which could be considered if one works to higher order 
in a n s (n > 3). 



Table 1: QCD input parameters. 
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Values 
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(4.5 ± 0.3) x 10~ 4 GeV 6 



134 36) 
average |2| 

(U ED 
HUM) 

EHUEilEML] 
□ 

01 EH 05) 



5. The LSR determination of the RGI fs in the static limit 

Analysis of the convergence of the PT series 

We study the effect of the truncation of the PT series on the 
value of fs from Eqs. (19i and (20 1. For a given value of 
p. = M b and a> c =3 GeV, we show the result of the analysis in 
Fig. [I] At the minimum (optimal) value, fs moves from 0.365 
(LO+NLO) to 0.414 (+N2LO) to 0.454 (+N3LO) GeV 3/2 , i.e. 
a change of about 13% from LO+NLO to N2LO and of about 
8.8% from N2LO to N3LO, which indicates a slow convergence 
of the PT series. We consider the N3LO contribution as a sys- 
tematic error from the truncation of the PT series. 




4 These values agree and improve previous sum rules results [4-7 40 41 1. 
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Figure 1 : a) r//-behaviour of for fi = and 0) c = 3 GeV for different truncations of 
the PT series . 



Analysis of the th and a> c stabilities 

We show in Fig. [2] the r w -behaviour of the result for a given 
value of p and for different values of cj c where the PT series 
is known to N2LO. A Testability is obtained for co c > 2 GeV, 



3 



while w c -stability is obtained for co c > 4 GeV. We consider as 
optimal and conservative values the ones obtained in the previ- 
ous range of <±> c values: 



f B (M b ) = (0.407 - 0.413) GeV 3/2 



(27) 



One should notice that contrary to some results obtained in the 
literature, we allow a larger range of <x> c -values from 2 GeV ( 
beginning of r//-stability) to 4 GeV ( beginning of w c -stability) 
in order to deduce a more conservative estimate. 




0.0 0.5 1.0 1.5 2.0 2.5 



Th [GeV ] 



Figure 2: a) r# -behaviour of for /i = M/ 7 and o» c - 3 GeV for different truncations of 
the PT series . 



Summary of the results for fy and error calculations 

At a given value of fi we estimate the errors induced by the 
QCD parameters compiled in Table [T] We summarize the re- 
sults of the analysis in Table [2] We show in Fig [3] the "sum 



Table 2: Central values and corresponding errors for f B in units of 10 3 xGeV 3/2 
from the LSR at different values of the subtraction point /j in units of GeV for 
Mb = 4804(50) MeV. The +(resp. -) sign means that the values of jb increase 
(resp. decrease) when the input increases (resp. decreases). The total error 
comes from a quadratic sum. 
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rules data points" at different values of the subtraction point /j. 
The error is large at small fi due to the bad behaviour of the PT 
series at low scale which confirms the scepticism of the authors 
of Ref. Il24l on the reliable extraction of f B at a such low scale. 
However, as we have shown in previous section, the conver- 
gence of the PT series improves obviously at larger scale which 
enables to extract f B with a reasonable accuracy. Fitting the 
previous data by an horizontal line or taking their average, we 
deduce the final value of the RGI universal coupling: 



f B = 0.405(25) GeV 
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Figure 3: fa at N2LO for different values of the subtraction point fi from Table m 
The green coloured region corresponds to the one where the error comes from a weighted 
average, while the non-coloured one delimited by the two green horizontal lines is the one 
where the error coming from the best determination has been taken. 



from which we can deduce the value of the static coupling eval- 
uated e.g. at M B =5.28 GeV: 



fs, a ,(M B ) = 0.587(2) ffi (37) /fi GeV 3/2 . 



The corresponding decay constant from Eq. ( 1 1 
f™ =228(l)„,(l%MeV, 



(29) 



(30) 



which is quite large compared to the value of fy = 206(7) MeV 
obtained from the full QCD theory [1| suggesting some large 
1 /Mb -corrections which we shall analyze in the next section. 



We consider the previous results in Eqs. (28 i to ( |30| > as im- 
provements of previous results in the literature 10314201 . 



6. 1/M b corrections and value of the decay constant f B 

Taking into account the mass-difference between the meson M B 

ex- 



and the on-shell quark mass Mb, the relation in Eq. (11 



pressed in terms of the RGI coupling in Eq. ( 12 1 becomes: 
C b (M b ) fl 



R b (M b ) M B 



(31) 



LSR expression of 1 /M b correction 5f^ 

The 1 /Mb corrections S fg to the HQET two-point correlator can 
be obtained by subtracting its expression in the full theory with 
the one of HQET in the limit M b — > oo. The l/M b correction to 
the physical decay constant reads (see e.g. EUll ): 



6ft = 
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where: 
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Up to order a s , it reads: 



lm6Yl PT (x) = - 
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8tt Ml 1 + x 
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where : x = a>/Mt and: 

F(x) = 2Li 2 (-x) + ln(x) ln(l + x) 
1 + x 

+ ln(l + jc) - 1 . 



1 + x 



ln(x) 



(35) 



The order N2LO a 2 PT correction to the spectral function can 
be numerically obtained by subtracting the complete expres- 
sion in full QCD obtained in 1 27 1 with the HQET asymptotic 



result in Eq. ( 18 1 and by using the relation in Eq. (32 1. In the 
same way, we estimate the N3LO PT corrections assuming a 
geometric growth of the PT series both in full QCD and HQET 
theories. 

The NP corrections read up to d=5 condensates |[T7ll20ll : 
$np(p) = 2a s (p.)(uu}(n) 



(36) 



M b 1 + cj I M b 



Analysis of the convergence of the PT series 
Like in the case of fg, we study the convergence of the PT se- 
ries. We notice that the a 2 s and a 3 s corrections are very small 
for th < 1 GeV, which can then be neglected. The analysis is 
shown in Fig. |4] 
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Figure 4: r^-behaviour of Sfg for fi = M/, and oj l = 3 GeV for different truncations of 
the PT series where the contributions of condensates up to d=5 have been included. 



Analysis of the t h and oj c stabilities 

We show in Fig. |5]the T^-behaviour of for different values 
of Ci) c and including the d — 5 condensates. We study the effects 
of the d — 5 condensates on the T//-stability for given two ex- 
tremal values of io c (beginning of t# and of w f -stabilites) . The 
analysis is shown in Fig. [6]from which we consider as optimal 
results the ones corresponding to the range: 

2.2) GeV" 2 , 



th - (1. 
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considering the fact that the inflexion point is not precisely lo- 
calized. 
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Figure 5: th -behaviour of <5/J for a = Mj } and for different values of uj c by including 
the d=5 condensate in the OPE. 



0.03 
0.02 
0.01 
0.00 
-0.01 
-0.02 



> 

a 

<N 
'<2 



-0.03. 





ai f [GeV] ii=M b 




-e-2, d=4 




-m-6, d=4 




- -A - 2, +d=5 








■ - ■ 6, +d=5 



0.0 0.5 1.0 1.5 2.0 2.5 
t h [GeV 1 ] 



3.0 



Figure 6: r//-behaviour of for fi = Mb and tL» t - 1.5 and 6 GeV for different trunc- 
tions of the OPE by the inclusion of the d=4 condensate or by the inclusion of d=4+5 
condensates. 



Results for 6 f^ 

The "sum rules data" of for different values of fi are shown 
Table [3] and in Fig. [7] where the main errors come from the 
localization of t# from 1.8 to 2.2 GeV 2 and the one induced 
from its corresponding cj c values. The errors from the QCD pa- 
rameters are negligible. Taking the average of different values, 



Table 3: Central values and corresponding errors for Sfl from the LSR at dif- 
ferent values of the subtraction point /i and for M/, = 4804(50) MeV. 



H [GeV] 


-6f B 10 j x[GeV 2 ] 


Error 10 3 x[GeV 2 ] 
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-0.27 


2.79 
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-1.19 


3.85 
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-0.66 


3.58 
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-0.40 


3.45 


M b 


-0.27 


3.45 
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-0.13 


3.32 
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3.45 



we deduce from the analysis the 1 /Mb corrections: 



Sfl = -0.48(1.37) x 1(T 3 GeV 2 . 



(38) 



One should notice that the mass correction from the central 
value of 5fg only decreases the value of fs by about 1 MeV 
which is negligible. 
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Figure 7: "Optimal sum rules data"' of 6f^ for different values /i. 



7. / B from HQET and from full QCD 



Combining the result in Eq. ( |38| with the one in Eq. ( [28} with 
the help of Eq. (31) , one obtains: 



jkqet _ 197 ( 13 ) MeV, 



(39) 



where we have added the errors quadratically. Notice that, un- 
like the full QCD case (TJ, we have not tried to extract an upper 
bound on f B from the positivity of the spectral function because 
of the undefinite sign of Sfg in Eq. (eq:fbphys). We can com- 
pare this result with the one obained in the static limit in Eq. 
( [30| ), where one can see that the main corrections are to the 
(M b /M B ) 3/2 ratio in Eq. 



(31 



The one due to 5 fg only de- 
creases the value of f B by about 1 MeV. One can also compare 
the HQET result with the one from the average of LSR and Mo- 
ment sum rules in full QCD JT|: 



(40) 



f q B cd = 206(7) MeV . 



Taking the average of the two sum rule predictions in Eqs. ( |39| l 
and ( |40| l , we deduce the final averaged prediction from QSSR: 

(41) 



jqssr = 2Q4(6) Me y 



8. Extraction of the b -quark mass from HQET 

One can extract the meson-quark mass-difference A using the 



ratio of the LSR obtained from Eq. ( 3 1 

(/>- T " AM ) 
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f2 e -T„AM 



AM = 



M\ -Ml 
M h 



(42) 



where Mb is the on-shell fo-quark mass. We notice that the 1 /M b 

corrections encoded in 5f\ which contributes about 0.9% of the 

value of /J, are also negligible for the extraction of AM though 

act to the th stability analysis of /J. We show the r//-behaviour 

of AM in Fig. [8] We show in Table |4] the different sources of 

errors on AM, where one can notice that the most important 

ones come from t c , the estimated a] and the mixed condensate 

cpjitributions. We show m Figi |9lthe /^-behaviour of different 
QSSR data points from which^we deduce the 



le average: 



Table 4: Central values and corresponding errors for AM in units of MeV from 
the LSR at different values of the subtraction point /i in units of GeV. The 
+(resp. -) sign means that the values of AM increase (resp. decrease) when 
the input increases (resp. decreases). The total error comes from a quadratic 
sum. 



A' 


AM 


tc 


a s 


J, 


M b 


(uu) 


<G 2 > 


Ml 


Total 


1 


941 


+ 124 


+ 10 


+21 


+ 1 


-7 


-5 


+22 


128 


2 


877 


+ 144 


+9 


+51 


+1 


-9 


-7 


+27 


156 


3 


841 


+ 147 


+9 


+58 


+1 


-11 


-8 


+27 


161 


4 


817 


+149 


+8 


+55 


+2 


-11 


-9 


+27 


162 


5 


798 


+149 


+8 


+64 


+2 


-12 


-10 


+27 


165 


6 


784 


+ 149 


+8 


+65 


+2 


-12 


-10 


+26 


166 



Using the previous value of AM, one can extract the on-shell 
fo-quark mass to order a 2 s : 



M 1 ^' =4871(30) MeV . 



(44) 



Using the known relation between the on-shell and running 
quark mass to order a 2 s (see e.g. loTl8l l32l ). we deduce: 



m b (mb) hqe ' = 4234(56) MeV 

2500, 

2000 



(45) 




oi, [GeV] ji=M b 

-a 1.5 

^3 
-- 4 

— 6 



0.5 1.0 1.5 2.0 
t h [GeV" 1 ] 

Figure 8: r//-behaviour of A for fi = Mb and for different values of c 




H [GeV] 

Figure 9: "Optimal sum rules data" of A for different values /i. 



9. m b {m b ) from HQET and from full QCD 

We can compare the previous value of the running mass with 
the one from heavy-light QCD spectral sum rules full QCD to 
order a] |T): 



AM = 853(63) MeV 



(43) 



m h 



(m b )i cd = 4236(69) MeV 



(46) 
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and from the T sum rules to order a] (2) : 

m b (m h ) r = 4177(11) MeV. (47) 
We can deduce the final averaged estimate from QSSR: 

m b {m b f ssr = 4181(11) MeV. (48) 

10. Summary and conclusions 

We have re-estimated fs and M b from HQET Laplace spectral 
sum rules to order a 2 s . Our results in Eqs. (39 1 and (45 1 are in 



good agreement with the ones from full QCD in Eqs. (40 1 and 



(46 1. We have combined different QSSR results including the 
value of m b (m b ) from the T systems in Eq. ( |47] i for deducing the 
final averaged estimate from QSSR in Eqs. (|41|l and (48 I. These 



results are in perfect agreement with some other determinations 
and with lattice results including nj — 2 or 3 dynamical quarks 
compiled in Table |3] EH56| . 



Table 5: Results for f'g and m~b(mi,) in units of MeV and comparison with lattice 
simulations using rif ■ = 2 1 53 54] and nf = 3 1551 .561 dynamical quarks, fp are 
normalized as f n = 130.4 MeV. 



Observables 



Methods 



Refs. 



h 

QSSR 

197(13) = 1.58(5)/, 

206(7)= 1.58(5)/* 

204(6) = 1.57(5)/, 

< 235.3(3.8) = 1.80(3)/, 

Lattice 

197(10) 

193(10) 

190(13) 

197(9) 

m b (m b ) 

QSSR 

4234(56) 

4236(69) 

4177(11) 

4181(11) 

Lattice 

4290(140) 



HQET 
full QCD 

Average 
full QCD 

ETMC 
ALPHA 
HPQCD 
FNAL 



S-meson - HQET 
5-meson - full QCD 
T - full QCD 

Average 

ETMC 



(this work) 
LD 

(this work) 

m 

[53| 
(51 
(55) 
[561 



(this work) 

LD 

in 

(this work) 



[53| 
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